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Abstract 

We study the thermodynamics of kinks in the $ 6 model using a Langevin code 
implemented on a massively parallel computer. This code can be used to study first 
order dynamical phase transitions which exhibit multiple length and time scales. The 
classical statistical mechanics of a 1 + 1-dimensional field theory reduces to a time- 
independent quantum problem in one dimension via the transfer integral method. 
Exact solutions of the Schrodinger equation exist for the $ 6 potential (unlike the 
case for $ 4 ) and can be used to check results from the simulations. The $ 6 model is 
also much richer than the <3> 4 model in terms of the variety of coherent structures and 
possible phases accompanying a phase transition. Specifically, we have calculated (in a 
one dimensional model) such quantities as the probability density function (PDF) and 
field-field correlation functions. These quantities help us understand the contribution 
to the specific heat from coherent structures such as domain walls (kinks) and other 
transformation structures as opposed to the contribution from lattice vibrations. We 
have calibrated our results against known exact solutions for limiting cases with very 
high accuracy. Having understood this problem, we are now extending our Langevin 
code to higher dimensions. 
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1 Introduction 



First order phase transitions are ubiquitous in nature, ranging from melting to struc- 
tural transformations in crystals. Due to the discontinuity in the order parameter 
at a first order transition, a simple symmetric double well potential is incapable of 
describing such a transition and an asymmetric double well or a symmetric triple well 
potential ( $ 6 model) is needed. A symmetric double well potential ($ 4 

model) is usually employed to describe continuous (second order) transitions. The 
$ 4 model and its attendant kink structure has been extensively studied in the litera- 
ture using techniques such as the path integral formalism , Langevin dynamics , 
etc. The <£> 6 model, appropriate for first order phase transitions, is much richer in 
terms of its kink structure ("textures" in the materials context) || @] || @ ■ However, 
unlike the $ 4 model, many aspects of the <3> 6 model, specifically the variation of the 
probability density function (PDF) of the order parameter j7|, correlation functions 
and structure factor within the context of Langevin dynamics, have not been studied 
(except in a specific case ||). 

Here we adopt a novel viewpoint to study the thermodynamics of a system based 
almost entirely on its PDF, i.e., if the PDF is known then the ground state eigen- 
function and eigenvalue are known, and therefore thermodynamic quantities, e.g., the 
specific heat, can be determined. (Correlation functions can also be determined but 
not directly from the PDF.) It is clear that the PDF must be determined to a very 
high accuracy. Our Langevin code, implemented on a massively parallel computer, 
computes the PDF to the required precision (in an appropriate temperature range). 
This can be easily checked for the $ 6 model because certain exact solutions of the 
Schrodinger equation with a $ 6 potential exist 00[|1O]- Indeed, we find spectacu- 
lar agreement between the computed and exact results for the $ 6 model. A direct 
calibration of the Langevin code in this manner against the <£> 4 model is not possible 
because the Schrodinger equation with a $ 4 potential has no exact solution. In ad- 
dition, the field configuration for a kink can be calibrated against the exact solution 
for the $ 6 kink. In this case too we find very good agreement between the Langevin 
code and the exact solution. 

2 The <p 6 Model 

The Landau free energy density for the <3> 6 model in one dimension is given by 

m ) 2 4 6 K ' 
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where to describe a first order (or discontinuous) phase transition it is necessary that 
B < and C > 0. In many physical cases (e.g., soft phonon mode driven structural 
transformations) A is chosen to be temperature dependent ||]||. This model also 
exhibits a second order (or continuous) transition for B > 0, C > at A = 0. To 
include domain walls (of nonvanishing width and energy) between various phases the 
Landau free energy is supplemented by adding the square of the field gradient (the 
Ginzburg term) 

F GL {$,&) = F L (§) + l -d& 2 , (2) 

where $' = d<&/dx and d > 0. The coefficient d also pertains to a physical quantity 
such as the soft shear modulus of a crystal || . For convenience we use dimensionless 
variables for the Ginzburg-Landau energy Fql, the field $, and the spatial variable 
x according to 

9B 3 

F G l = AF GL , A = - , (3) 

/ 3£\ 1/2 

$ = 7$ , 7 = (-—J , (4) 

Thus, in dimensionless form we get 

F GL = ir$ 2 - $ 4 + $ 6 + $' 2 , (6) 

where the prime denotes differentiation with respect to x, and 
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can be interpreted as a dimensionless temperature. The <J> 6 potential (|1|) has three 
minima that occur at 

5 = 0, $ min = ±$ (l + ^) 1/2 , (8) 
*o = ^7, ,= (l-jr) , (9) 

and two maxima that occur at 

$ max = ±$ (1 - ^) 1/2 • (10) 
The value of the potential (Landau free energy) at the extrema is given by 



F ± = F (lTzf(l±2z) , 
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(11) 



where F + represents a free energy maximum, F_, a free energy minimum, and F = 
(2/27) A. For the dimensionless free energy Fql and field $ one just sets 7 = A = v — 1 
in Eqs. (§£!]). 

Physically the central minimum represents the high temperature (parent) phase 
while the two side minima represent the two variants of a fully developed low tem- 
perature (product) phase. Similarly, the two maxima refer to the two variants of a 
partially developed product phase. For r < there are only two minima and the free 
energy behaves as an effective $ 4 model. For < r < 4/3 there are three minima. 
At r = 1 there are three degenerate minima which correspond to the first order tran- 
sition point, i.e., the parent and product phases coexist in equilibrium. For r > 4/3 
there is a single well corresponding to the fact that only the parent phase is stable. 
At r = 4/3 there is only one minimum at $ = and two points of inflection at 
±(-5/2C) 1/2 . 



3 Exact Kink and Domain Wall Solutions 

The equilibrium field configuration is determined by minimizing the total free energy. 
The Euler-Lagrange equations after two integrations lead to 

with cf) = $ 2 and the boundary conditions F = lim^ioo Fl( x ), (Fl{x) > F ), 
lim^-i-oo $' = 0. The following four kink and domain wall solutions exist. 

1. For t < 1 and Fq < a kink solution between the two product variants is 
given by 

$ min asinh/3x 

(X) = ~o o 2 \V2 ' (13) 

r + a 2 sinh (3x ) 



where 

, (3 = $ mm (3$Ln " l) • (14) 



a = $ T 



2$L„ - 1 



This kink (Fig. 1) corresponds to the $ 4 kink for r < 0. 



2. For 1 < r < 4/3 and F > a pulse soliton between the parent and either product 
variant (with the parent phase in the middle) is given by 



Figure 1: Kink solutions of type 1 for r = .99, .5, — .5 (asymptotically, the bottom, 
middle, and top curves). Note the "sticky" behavior of the solution near = for r = .99. 



This pulse solution is depicted in Fig. 2. The total free energy for solutions 1 and 2 
is given by 

Ftot = \U ~ (r - 1) In ( ^ + *" iw ) ) • (16) 

3. For < t < 1 and F = a pulse soliton between either product variant and the 
parent phase (with either product variant in the middle) is given by 

Q(X) = — 77T , (17) 

/ / <T>2\ o \ 1/2 ' v ' 

(l+(l-||)sinh 2 ^/2) 

where $23 = (1/2) (1 =F (1 — t) 1 ^ 2 ). This solution is shown in Fig. 3. The total free 
energy for this pulse solution is given by 

4. For r = 1 and Fq = a "half" kink between either product variant and the parent 
phase (at the transition point) is given by (Fig. 4) 

$(*) = i=(l + e-*)" 1/2 , (19) 

with total energy F tot = 1/8. This is the limiting case of 3 (matching with half a 
pulse of type 3). Note that the above four solutions are known in the literature in a 
different form (and materials context) ||[Q||fll][l^]. 
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Figure 2: Pulse solutions of type 2 for r = 1.0001, 1.1, 1.3. (asymptotically, the bottom, 
middle, and top curves) 




Figure 3: Pulse solutions of type 3 for r = .1, .5, .999 (the top, middle, and bottom peaks 
respectively). 
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Figure 4: Half kink solution of type 4. 



4 Langevin Dynamics on a Parallel Machine 

Solution of field theoretic Langevin equations is particularly convenient on a mas- 
sively parallel computer because little inter-processor communication is involved and 
because the large memory available enables the use of large lattices. For the results 
reported here, lattice sizes ranged from several thousands to hundreds of thousands of 
lattice points. The system size was kept much larger than the field correlation length 
(at least a factor of ten, but much larger typically). Periodic boundary conditions 
were used for convenience. 

The Langevin equation for the $ 6 model is 

4$ = %J} - 77<9 t $ + $(1 - $ 2 ) + $ 5 + F(x, t) , (20) 

where the viscosity rj and the Gaussian white noise F are related by the fluctuation- 
dissipation theorem: 

F(x, t)F(x', t')) = 2t 1 i3- 1 5(x - x')5(t - t') . (21) 



The lattice versions of the above continuous equations were then solved using standard 
techniques ]13| . Random initial conditions were driven to equilibrium and the results 



sampled in time thereafter to yield time averaged PDF's, etc. The use of the Langevin 
technique for obtaining thermodynamic quantities is straightforward and remarkably 
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accurate. Moreover, structures such as the various kink solutions can be clearly 
identified and information about real time dynamical quantities such as the temporal 
correlation functions is also available. 



5 Results 



5.1 Probability Distribution Function 



The exact solution for the PDF is given by (the ground state wave function [|9| [ 10 1 
squared for the quantum mechanical problem is the PDF) 



#jj = Nexp { -(3 



1 ( dC \ 1/2 ^ , B 



(22) 



where N is a normalization constant. The corresponding ground state energy [|9[ [|10 

In Fig. 5 we compare the probability density (PDF) function computed from a 
Langevin simulation (diamonds) with the exact solution (dashed line). The choice of 
parameters is A = B = d = 1, C = 0.17, and 1/(3 = 0.144146. As is apparent, the 
agreement is excellent, and is indicative of the very high accuracy of the Langevin 
simulations. The results of our computations are, for practical purposes, numerically 
exact. Note that such a comparison is not directly possible in the case of the $ 4 
model since there are no exact solutions in that case. 

We have studied the variation of the PDF with temperature for several temper- 
atures ranging from well below the "transition" temperature to well above it. In 
general, the PDF exhibits a three peak structure. However, at a certain temperature 
the three peaks have the same height. Above this temperature the PDF is charac- 
terized by a dominant central peak whereas below this temperature there are two 
dominant side peaks. In higher dimensions, this is characteristic of a first order phase 
transition. Since the PDF contains all essential thermodynamic information it is very 
important to be able to compute it accurately: our approach provides a simple and 
accurate method for computing thermodynamic quantities such as the specific heat, 



internal energy, etc. A detailed description is now in preparation [14 



5.2 Field Configurations 

A low temperature field configuration is shown below in Fig. 6. The kinks are few 
in number but well defined. At higher temperatures the number of kinks increases 
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Figure 5: Comparison of the exact (dashed line) and numerical (diamonds) PDF. The 
agreement is excellent (parameter values are given in the text). 




2000 4000 6000 8000 10000 12000 14000 16000 

Figure 6: A sample field configuration at 1/(3 = .2 with all other parameters the same as 
in Fig. 5. 
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but their shape is smeared by thermal noise. At still higher temperatures it becomes 
impossible to distinguish these nonlinear structures from nonlinear phonons. At mod- 
est temperatures, the shape of the kink computed from Langevin simulations can be 
compared with the exact solutions described in Sec. 3 (for example, the kinks in Fig. 
6 correspond to the kinks of Fig. 1 with r ~ .91). Here we simply report that the 
results from the simulations match well with theoretical expectations. Details will be 
given elsewhere ||14| . To summarize, we have two different checks, namely the PDF 
and kink shape, which both show very good accuracy and provide a high level of 
confidence in our simulations. 



5.3 Correlation Functions 

The absolute value of the location of the side minima is the order parameter for 
the first order phase transition. Interestingly, the topological charge associated with 
a kink interpolating between these two minima at the phase transition point (T c , 
three degenerate minima) turns out to be precisely equal to the order parameter 
0. In addition to the value of the order parameter, its spatial correlations as well 
as correlations of its intensity are often of interest in studying a phase transition 
IJ0. The correlation functions are particularly interesting because they describe 
the behavior of systems which are nearly ordered but do not undergo sharp phase 
transitions at any finite temperature. Following the path integral (transfer operator) 
procedure 00, the correlation functions are expressed in terms of the eigenvalues 
and eigenstates of the transfer operator as follows: 



C x {x) = £|(^ n |$|^ )| 2 exp 

n 

C 2 {x) = £|(^„,H$| 2 |^o> 



-/3-(e n - e ) 
a 



2 

exp 



-f3-(e n - e ) 
a 



(24) 
(25) 



where p = l/k B T, S |$| 2 = |$(x)| 2 - (|$(x)| 2 ), and a is the lattice constant. 

For x ^> £, the lowest excited state coupled by the matrix element determines 
the behavior of the correlation functions. In other words, the eigenvalues set inverse 
correlation lengths. For T ~ T c , near degeneracy in eigenstates is reflected in an 
increased range of correlation (tunneling) 00. At large distances C\(x) and (^(x) 
are dominated by the state with smallest eigenvalue for which the corresponding 
matrix elements are nonvanishing (excluding the n = state). The correlation lengths 
for C\ and C 2 are, respectively, 

y * -(ci-eo), (26) 
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Figure 7: The (unnormalised) correlation function C\ at four different temperatures plotted 
on a logarithmic scale. The parameters are the same as in Fig. 5, with only the temperature 
being varied: from top to bottom = .3, .35, .4, .45. 
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(e2 - e ) 



(27) 



£i is proportional to the average separation between neighboring kinks, which is 
the distance over which the field remains correlated. £i grows exponentially with 
decreasing temperature, and £i^ooasT^0, when no kinks remain in the system. 
The energy density correlations are usually short range. The static structure factor 
(or the equal time correlation function) is the Fourier transform of C\, and is given 
by 

S(q) = — [ dxe iqx (<£>(0)$(a;)> . (28) 

The field-field (or order-parameter-order-parameter) correlation function C\ for 
four temperatures is shown in Fig. 7. The exponential decay is apparent. The 
correlation length is given directly by the slope of the correlation function plotted 
on a logarithmic scale, while the average domain size in the system is obtained from 
the first zero crossing. The corresponding structure factors (fast Fourier transform of 
C\) are also easy to compute but we do not display them here. More details on the 
correlation functions, their exact and semi-exact calculation, and comparison with 



numerical results will be given elsewhere [14 
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6 Conclusion 

In conclusion we restate some key points. First, the <J> 6 theory has sufficient structure 
to describe first order phase transitions especially relevant in the materials context 
(shape memory alloys) A consequence of this complexity is the appearance 

of several coherent nonlinear structures. The thermodynamics of the theory can be 
profitably studied via both the transfer operator method and Langevin simulations. 
The remarkable occurrence of some exact solutions in the analytic transfer operator 
approach for the $ 6 theory allows B for a strong check on the simulations. The 
very accurate determination of the PDF via our simulations implies that this maybe 
a convenient window for a study of the thermodynamics of such systems. Finally, 
extension of the Langevin method to higher dimensions, and other classes of quasi- 
exactly solvable potentials, is simple and we expect to present our results for both 
the two and three dimensional cases soon. 
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